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1 Introduction 

The quest for a better understanding of the properties of nuclear matter under ex¬ 
treme conditions (such as those created in relativistic heavy ion collisions) has led 
to a number of theoretical challenges. One of them is to explain the success of the 
hydrodynamic description of quark-gluon plasma (QGP) evolution at time scales of 
order r ~ 0.5 — 1 fm/c [I1I2]- Immediately following the nuclear collision the system is 
very far from equilibrium, yet rapidly reaches a state where hydrodynamics appears 
to apply. This process is fast in the sense that it occurs on a time scale less or equal 
the inverse local temperature, i.e. thTh < 1. 

The fact that the hydrodynamic description works well with a very low value of the 
shear viscosity to entropy ratio (rj/s) and vanishing bulk viscosity suggests that the 
plasma state may be viewed effectively as a strongly coupled, conformal fluid. This has 
motivated studies aimed at understanding how such a hydrodynamic description could 
emerge in a case which is amenable to theoretical studies - the A/” = 4 supersymmetric 
Yang-Mills theory (SYM). In this example one has the option to apply gauge/gravity 
duality [3] to model the approach to equilibrium. Such studies were initiated in [1], 
where the case of Bjorken flow (discussed in more detail below) was considered. Janik 
and Peschanski showed in that context that the hydrodynamic description does indeed 
emerge at late proper times. This work was followed by numerous articles devoted 
to the description of hydrodynamic states in the context of gauge/gravity duality. In 
particular, second order transport coefficients were computed [SI E], and many new 
insights into the meaning of relativistic hydrodynamics were gained [HEiEiiiniin]. 

To describe far from equilibrium states of A/" = 4 SYM it is necessary to resort 
to numerical calculations. A decisive step opening this field of research was made 
by Chesler and Yaffe [12] who devised a very effective numerical scheme for solving 
Einstein equations in asymptotically AdS spaces based on characteristic evolution. 
This was soon applied to the case of Bjorken flow [T3|, which is a particularly attrac¬ 
tive setting, since it was an important model used to understand basic features of 
QGP evolution (such as entropy production), and at the same time is simple enough 
to implement easily in the context of gauge/gravity duality. The results described 
here [H] follow from a modification of this scheme. A different method was used 
in [T31 [TB] . where some basic features of the approach to hydrodynamics were stud¬ 
ied. Most importantly, it was found there that the system reaches the hydrodynamic 
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regime quickly (in the sense described earlier). Another physically important conclu¬ 
sion from these papers (earlier noted in [13]) was that hydrodynamics works very well 
already at a time when pressure gradients are large. Thus the process of reaching this 
stage of evolution is often referred to as hydrodynamization instead of thermalization. 

The numerical studies Ha HE] used a different numerical scheme from Ha HE], 
which was however limited by the fact that only a few (29) consistent initial states 
were known. The study jH] reported here adapted the approach of [13] in a way which 
allows an arbitrary number of initial conditions to be analysed, making it possible to 
look for generic features. In our work we looked at 600 initial conditions randomly 
generated on the gravity side of the duality. In the context of gauge/gravity duality 
there is a natural characteristic of the initial state called the initial entropy [15] (to 
be dehned precisely below). Observables such as hydrodynamization time depend in 
particular on this quantity. 

In this note we will focus on two natural questions 

• Is hydrodynamization generically a fast process? 

• Are there any universal physical characteristics of hydrodynamization? 

The answer to the hrst question seems to be positive and conhrms previous investiga¬ 
tions US] HE]; regardless of the values of the initial entropy chosen the system reaches 
hydrodynamic description on time scales of the order of the inverse local tempera¬ 
ture. The second question is more subtle. In the sample of initial data analysed in our 
study the hydrodynamization time appears not to be correlated with initial entropy 
(in contrast to [TE], where such a correlation seemed to be present on the basis of a 
smaller sample of initial states). However, when looking at the energy density at the 
time when hydrodynamic evolution starts, there appears to be a linear correlation 
with initial entropy. 

One should also mention the work [13 HE] where the process of isotropization was 
considered in a similar spirit. Recently a technically different (closer to the original 
formulation of the problem CZI) but conceptually similar studies appeared in [T9] . 

2 Supersymmetric plasma and gravity 

As discussed in the introduction, TV" = 4 SYM is amenable to quantitative studies in 
the strongly coupled limit. In the absence of methods which could be used in QCD 
for the study of strongly coupled, real time dynamics, this theory has become a fo¬ 
cus of much attention. This theory shares some features of QCD (especially at high 
temperature), but is rather different from it. It contains, apart from gluons, 6 real 
massless scalars and 4 Majorana massless fermions, all in the adjoint representation 
of the U{N) gauge group. The theory is known to be conformal even at the quantum 
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level. As mentioned in the introduction, we work within the AdS/CFT correspon¬ 
dence [3] which becomes an effective computational tool in the’t Hooft limit N ^ oo 
and A = OO) where quantum and stringy corrections on the gravity side can 

be neglected. 

There are important similarities and differences between SYM and QCD which one 
has to keep in mind. Among the similarities are the existence of the deconfined phase. 
Also, in the perturbative regime at T > 0, both theories have been shown to behave 
similarly, with the difference coming mostly from the different number of degrees 
of freedom |2D]. The most crucial difference is that A/” = 4 SYM has a vanishing 
beta function, which implies that it is not confining, has no finite temperature phase 
transition and has an exactly conformal equations of state. 

The modeling of nuclear collisions is an extremely complex task. To reduce the 
complexity of the problem we adopt strong symmetry assumptions introduced by 
Bjorken for the description of matter following a heavy ion collision. The dy¬ 
namics of the system is assumed to be independent of boosts along the longitudinal 
(collision) axis and is independent of transverse coordinates. In proper time-rapidity 
coordinates t = t cosh y , z = t sinh y this reduces to the statement that observables 
dependent only on proper time r. This approximation becomes exact in the limit of 
an inhnite energy collision of inhnitely large nuclei. 

Two physical quantities of interest to us will be the energy momentum tensor and 
entropy. In the present circumstances first one takes the form 


= Diag{e(r),pL(r),pr(r),PT(T)} , 


( 1 ) 


where e(r) = Pl{t) -|- 2pt{t) as required by conformal symmetry. The energy density 
defines local effective temperature by the relation 

£(r) = ?/VVT(r)‘‘ . (2) 

This is the temperature of an equilibrium system with the same energy density. It 
can be shown that for late times the dynamics is governed by the equations of hydro¬ 
dynamics; up to third order the effective temperature follows 


T(r) 


A n _ 1 , -l + log2 

Gtt (Ar)^'^^ SGtt^ (Ar)"^'^^ 
-21 + 271^ + 51 log 2-24 log^ 2 > 
19447r3 (At)^ ^ 


+ 


(3) 


The energy scale A appearing in Eq. (|^ depends on the initial conditions chosen 
and it is the only trace of initial state information contained in the hydrodynamic 
expansion. It also sets the scale for the energy density at hydrodynamization. 
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The calculation strategy follows usual lore of holography PI- States in the bound¬ 
ary theory correspond to asymptotically AdS geometries in the bulk. For example 
an equilibrium, hnite temperature, deconhned plasma state corresponds to a static 
(planar) black hole in the bulk. The Hawking temperature of this black object is 
interpreted as the temperature in the dual A/” = 4 SYM. Extending this notion to 
the out-of-equilibrium situation we assume that non-equilibrium plasma states corre¬ 
spond to geometries with non-static horizons. Such geometries are assumed to poses 
an event horizon, but there are reasons to believe ^31I7F] that the physical notions 
(such as entropy for instance) should be associated with apparent horizons. 

This conjecture allows us to extend the notion of entropy to non-equilibrium states 
by the Bekenstein-Hawking relation, which with our normalization translates to 

g _ OaH 
71 

where oah is apparent horizon area [23]. By the area law theorems this quantity is 
non-decreasing and agrees with thermodynamic dehnition for late times. 

In calculations performed here we took 600 different initial states described by 
randomly generated initial geometries - with each of these we associate an initial 
entropy as dehned above. We then evolved these geometries according to Einstein 
equations well into the hydrodynamic regime. From the rules of the holographic 
correspondence we are able to read of the relevant physical observables i.e. energy 
density e(r) and entropy S{t). An important assumption on the initial conditions 
is that e(0) ^ 0, which allows us to normalize the initial effective temperature as 
T(0) = I/tt. For more technical details of construction of initial geometries and 
solving for the time evolution we refer to the original paper |14j . 

3 Results and conclusions 

In order to present quantitative results on the hydrodynamization process we need 
to give it a precise dehnition. The approach to hydrodynamics can be observed by 
monitoring the pressure anisotropy 


^ ^ Pt-Pl 
e 

It is convenient to measure time in units of inverse local temperature, that is, to 
use w = tT{t) as a parameter of evolution. The pressure anisotropy can then be 
expressed as A(tc) = Qf{w) — 4 in terms of the function (hrst introduced in |15j j 


/H 


r dw 
w dr 


( 6 ) 
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which for dimensional reasons is independent of A. At large times (large w) this 
function attains a universal (hydrodynamic) form determined by an infinite 

set of transport coefficients, which up to the third order reads 


fniw) 


2 1 l-log2 15-27r2-451og2 + 241og^2 

3 971W 277r‘^w^ 9727r^t(;^ 


(7) 


The beginning of the hydrodynamic stage might now be defined as the value of w 
(hence proper time r) when the difference between the actual f{w) and the hydrody¬ 
namic form fniw) is less than some arbitrary small number; for example 


/g(w) 

fH 


< 0.05. 


( 8 ) 


While this definition involves some arbitrary choice, varying this criterion within rea¬ 
son leads to no appreciable change in the calculated value of the hydrodynamization 
time. 



Wth 


Figure 1: Histogram of hydrodynamization times Wth in units of effective hydrody¬ 
namization temperature. 

Figure shows a histogram of hydrodynamization times obtained this way. Re¬ 
gardless the initial entropy, the hydrodynamization time is of the order if the inverse 
hydrodynamization temperature. On the average Wav = 0.57. It is instructive to com¬ 
pare it to the estimation from the RHIC data; T = 500 MeV and r = 0.25 fm/c gives 
■inRHic = 0.63 which is very close to theoretical prediction. The pressure anisotropy 
at hydrodynamization is found to be quite high: A 0.35. 

The other quantity we focus on here is the energy scale A which sets the scale for 
the hydrodynamic cooling. This quantity is obtained by fitting the tail of the data 
to Eq. 
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Figure 2: Hydrodynamic energy scale A in units of initial effective temperature. 


Results for this scale are shown in figure For an intermediate range of entropies 
below S ~ 0.3 a strong linear correlation is observed. For initial entropies larger than 
0.3 the correlation is lost and “chaotic” behaviour develops. 

In conclusion, through our analysis of the time evolution of a bulk sample of 
different non-equilibrium initial states we find signihcant evidence to claim that the 
early hydrodynamiztion of [T3] is a generic process. Our analysis supports previ¬ 
ous hndings that hydrodynamiztion is different from thermalization in the sense of 
non-negligible pressure gradients being present and well described by hydro. It also 
suggests that entropy production during the hydrodynamic stage of evolution is not 
negligible despite the low value of rj/s. 

Our simulations suggest that, at least in some regions of initial state parameters, 
there might exist characteristic regularities, reflecting the nature of collective non¬ 
hydrodynamic degrees of freedom [2H]. At this point an important question to what 
extent is this correlation a consequence of the strong symmetry assumptions imposed 
and to what extent it reflects the true nature of the process. 
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